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Due to a transient quantum interference during a 
wavepacket collision with a potential barrier, a particular mo- 
mentum, that depends on the potential parameters but is 
close to the initial average momentum, becomes suppressed. 
The hole left pushes the momentum distribution outwards 
leading to a significant constructive enhancement of lower 
and higher momenta. This is explained in the momentum 
complex-plane language in terms of a saddle point and two 
contiguous "structural" poles, which are not associated with 
resonances but with incident and transmitted components of 
the wavefunction. 



The traditional formulation of quantum scattering the- 
ory in terms of an "S-matrix" assumes that only the re- 
sults of the collision can be observed at "asymptotic" 
distances and times, but that the collision itself cannot 
be observed. This perspective is justified to analyze the 
products of standard collisions of atomic or molecular 
beams. But the iS-matrix approach is not enough to de- 
scribe modern experiments where the collision complex 
can be observed by means of femtosecond laser pulses or 
"spectroscopy of the transition state" fl. Also, in quan- 
tum kinetic theory of gases accurate treatments must 
abandon the "completed collision" approximation and 
use a full description, e.g. in terms of Moller wave oper- 
ators as in the Waldmann-Snider equation and its gen- 
eralizations for moderately dense gases In any case, 
it is important to understand the collision process itself 
to control or modify the products. This has motivated a 
recent trend of theoretical and experimental work to in- 
vestigate the details of the collision itself, and not only its 
asymptotics. In particular, a quantum effect has been re- 
cently described by Brouard and Muga ||f| in which the 
probability to find the particle with a momentum above 
a given value is larger, in the midst of the collision, than 
the quantity allowed classically by energy conservation. 
The effect belongs to a group where the conservation of 
classical energy seems to be violated. Well known exam- 
ples are the tunnel effect, or in general the non vanishing 
probability to find the particle in evanescent regions be- 
yond the classical turning points. 

This transient enhancement of the momentum tail may 
in principle be observed by collisions of ultracold atoms 
with a laser field that can be turned off suddenly in the 
time scale of the atomic motion ||, and implies as a 
macroscopic consequence deviations from the Maxwellian 



velocity distribution ||] . We initiated the research of the 
present work looking for conditions that increase the ef- 
fect and favor its observability. In so doing we have found 
an unexpected regime where the effect is much higher 
than in previously studied cases. In this letter we shall 
describe such a regime and analyze its physical origin, 
namely a transient interference between transmission and 
incidence components of the wavepacket. Let us first re- 
view briefly the main aspects of the classically forbiden 
increase of high-momenta. Brouard and Muga have stud- 
ied several examples where the quantity 



{\i>{p',t)\ 2 -\^(p',Q)\ 2 }dp' 



takes on positive values for positive potentials (the cor- 
responding classical quantity is negative or zero due to 
energy conservation) |||]]. An important aspect of this 
effect is its transient character, G q < before and after 
the collision. The effect is also generic because the 
stationary components of the wavepacket have, in mo- 
mentum representation, a tail due to the resolvent which 
is always present in the Lippmann-Schwinger equation. 
This tail goes beyond the maximum value allowed by the 
conservation of energy. 

For a Gaussian wavepacket colliding with an infinite 
wall, maximum values of G'^ lax (p,t) ~ 0.05 have been 
reported Also a "delta" potential was used Q to 
analyze the influence of the opacity of the barrier. For 
the cases examined, an increase of G q with the opac- 
ity was observed up to a saturation level where the in- 
finite wall results were recovered [Q. This suggested 
that the observability of the effect would improve with 
strongly opaque conditions. In a complementary study 
we have systematically varied the spatial variance of the 
wavepacket, 8 X , and the height of a square barrier, Vq, 
for a fixed average initial momentum p c . We have found, 
contrary to previous expectations, that the maximum ef- 
fect corresponds to energies well above the barrier and to 
large values of 5 X . In this regime the barrier is not at all 
opaque and essentially the full wave is finally transmit- 
ted. Moreover, G£, oa . is as large as 0.27. 

The numerical effort to perform these calculations by 
ordinary propagation methods (such as the split operator 
method) is rather heavy, since large values of S x and the 
need to discern fine details of the momentum distribution 
require an extense and dense grid. In fact for very large 
values of 8 X this numerical route has to be eventually 
abandoned. But even if one gets numerical results with 
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enough computer power, they will not provide any expla- 
nation of the unexpectedly high G q values. Fortunately 
these two difficulties can be overcome with an approx- 
imate analytical solution. Here we shall sketch its ob- 
tention, a more detailed account will be given elsewhere. 
First the momentum representation of the wavefunction 
is expressed using the basis of stationary eigenstates of 
H, \p + ), corresponding to incident momentum p' , and 
energy E' — p' 2 / (2m), 
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If the initial state at time t = does not overlap with the 
potential and has negligible negative momentum compo- 
nents we can write 
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To facilitate the treatment of the integral in the p'- 
complex plane we may now extend the lower limit to 
—00 using the analytical continuation of (p\p + ), p' > 0, 
over p' < (and later over the whole complex plane) . 

For a square barrier of height Vb and width d, centered 
at the coordinate origin, the delta-normalized stationary 
wavefunction with incident momentum p' is 
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x < -d/2 
-d/2 < x < d/2 
x > d/2, 

(3) 



where I = I, k' = p 1 /% and k" = yjp' 2 - 2mV /h. The 
cocfficientes R,C, D and T are determined by continu- 
ity of the wavefunction and its derivative. The momen- 
tum representation (p\p' + ) will correspondingly have five 
terms. The terms with /, R and T have poles in the 
p'-complex momentum plane at 



Pi 



p + iO 
p' R = -p- i0 

Pt=P- i0 ! 



(4) 



while the terms with C and D do not have these struc- 
tural poles H , which are not related to resonances or to 
the potential profile. The four functions R, C, D and T 
present an infinite series of resonance and anti-resonance 
poles in the third and fourth quadrants due to the zeros 
of a common denominator 

Q(p') = cos (k"d) -t^ + ^j sin (k"d) . (5) 

(In particular T(p') — exp(— ik'd)/Vl(p').) The condi- 
tions examined in this work correspond however to "di- 
rect scattering" , where these resonance singularities do 
not play any significant role. 



The initial state is taken as a minimum-uncertainty- 
product Gaussian centered at the position — aS x , a > 0, 
with average momentum p c , 
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This expression and the momentum representation of (J3|) 
are inserted in (||) to obtain five integrals. The full 
treatment of the resulting integrals follows closely ref. 
0. The integrals with C and D may be evaluated with 
the steepest descent method for large values of 8 X . The 
steepest descent path (SDP) is a straight line with slope 
—th/(2mS x ), with a saddle point close to p c in the midst 
of the collision. We shall always assume that the slope is 
small so that when the integration contour is deformed 
along this path it "cuts" the resonance poles of the fourth 
quadrant far from the real axis, and their residues can be 
neglected ("direct scattering" conditions). 

Because of the interference between the saddle and the 
structural poles, the simple steepest descent treatment 
valid for C and D is not valid for the other terms. A 
uniform expression for a smooth treatment of the pole 
crossing of the SDP is provided by the w-function, w(z) = 
e~ z erfc(-iz), which may also be defined by its integral 
expression H 



w(z) = — / du — 
in J r _ u - 



(7) 



where T_ goes from — 00 to 00 passing below the pole. 
Since we are interested in wavepackets with energy well 
above the barrier maximum the "reflection term" with R 
may be neglected. The remaining contribution is 



ipiT = ih 1 / 2 hr 
where 
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and 
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The functions gi(j>') and grip') present structural poles 
at p'j and p' T ; in addition grip') has resonance and anti- 
resonance poles. 

The SDP is a straight line with the same negative slope 
as before, passing through the saddle point, 
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To integrate (g), the contour is deformed to the SDP 
passing over the saddle. The same reasons to neglect the 
residues from the resonance poles in the C and D terms 
are now applicable. To introduce the w-functions, the 
integrand must be put in the form (0) . We complete the 
square in (|)J) and use the change of variable 
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to obtain 
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where g{u) = g(p'(u)), and 
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We may retain the main contribution from gi from 
its behaviour near the pole by approximating giiu) sa 
1Zi/(u — uj), where IZi is residue of gi{u) at the point 
u = uj = (p'j — /)//, and similarly for gp. For an approx- 
imate expression of {p\ip), and considering that the wave 
is much more extended in space than the barrier we may 
neglect the contribution from C and D and retain only 
the incidence and transmission terms, 

x \w{ui) + T(p)w(-u T )} = i^fa t). (12) 

A more precise expression including a term i/jrcd arid 
corrections to the zeroth order ip® T is given elsewhere 
and allows to obtain the wave function and G1 nax ac- 
curately for large values of S x with small computational 
effort. However ( |l2|) captures the essential, and provides 
a simple, explanatory picture of the phenomenon we want 
to discuss. 

Fig. 1 shows the distribution of momenta | (p\ip(t)) | 2 
for different instants of time, from the initial one to a time 
after the collision has been completed, passing through 
the instant where G q = 0.27 is maximum. In all figures 
the numerical values correspond to collisions of ultracold 
Rubidium atoms with an effective laser potential. The 
observed behaviour does not have a classical explana- 
tion. Note that the wavepacket is considerably broader 
than the barrier. A classical ensemble of particles with 
the same Gaussian phase-space (Wigner) distribution as 
(||) would only be slightly deformed due to the small frac- 
tion of particles located on the barrier at a given time, 
and would keep the maximum at the average momentum 



p c . Moreover, there could not be any spectacular accel- 
eration or deceleration as the one seen in the two peaks 
of the quantum distribution. We shall see that the zero 
of the quantum momentum distribution, which forbids 
in this case the initially dominant momentum p c , is due 
to a destructive interference whereas the two new peaks 
correspond to momentum regions of constructive inter- 
ference. 

In Fig. 2 the Argand diagrams of the two terms are 
represented, namely the imaginary versus the real parts 
obtained by varying p at equal intervals. Each lobule 
corresponds to one of the terms. The "motion" as p in- 
creases begins close to the origin, downwards in both 
diagrams. The left peak of the momentum distribution 
corresponds to the zone where the two moduli increase 
together and are approximately in phase. After the de- 
scending motion there is a fast, aproximately circular mo- 
tion where the phases become opposed (destructive inter- 
ference). Finally, the two curves meet again in phase in 
the upper part of the lobules, this momentum interval 
corresponds to the right peak of the momentum distri- 
bution. The described behaviour is essentially due to the 
two w-functions w(ui) y —w(—ut), as shown in Fig. 3, 
where the two Argand diagrams of the two w-functions 
and of the factors that multiply them are represented be- 
tween the momenta of the two maxima. Clearly the efect 
of the factors, whose phases remain essentially constant 
around 7r, is simply to invert the two lobules of the ws. 
The fast motion of the w-functions is due to the pass of 
the two contiguous structural poles uj = (p + i0)/f and 
u t = (p — *0)// near the saddle point at u = 0. A sweep 
from smaller to larger values of p moves the couple of 
poles along the real p' axis from left to right, while, for 
fixed t, the saddle point and the steepest descent path 
do not depend of p. Since ui ~ ut we can write, using 
the relation between w-functions of argument of opposite 
sign, see (|), 



w{ui) 



w{—ut). 



(13) 



During the collision, the saddle point is very close to the 
real axis of the p'-plane, only slightly below in Fig. 4, 
and the slope of the SDP is very small. This means that 
the difference between the two ^-contributions is essen- 
tially a real exponential, which implies a "simultaneous 
motion" , with equal maginary parts, along the two lob- 
ules of the Argand diagram. The rapid variation of the 
phases of the ws when passing close to the saddle point 
follows from its integral expression (Q). When uj passes 
close to u — and close to the real axis of the it-plane, 
the denominator is essentially real and changes its sign 
quickly, so there is a rapid change by 7r in the phases 
of w(ui) and -to(-ur). The phase oposition alone does 
not explain however why the interference is totally de- 
structive. It is also necessary to have equal moduli of the 
two incidence and transmission terms of (|^|) for an ex- 
act cancellation. Actually the equality is obtained only 



3 



transitorily, since before and after the collision only one 
lobule remains, the one for incidence before the collision, 
and the one for transmission after the collision. Along the 
collision the incident lobule decreases and the transmis- 
sion one grows, until they equilibrate and give a perfect 
cancellation and the two constructive interference zones 
of Fig. 1. 

By changing the barrier height the fases of the factors 
that multiply the ws change, the lobules rotate with re- 
spect to each other, and one of the two in-phase regions 
grows while the other diminishes, so that the two peaks of 
the momentum distribution become asymetric, see Fig- 
ures 5 and 6, where the momentum distributions and the 
corresponding lobules of the Argand diagrams are shown, 
compare also with Fig. 2. Note that these factors do 
not depend on time and therefore the angle between the 
lobules remains constant throughout the collision. This 
means that the positions of the maxima and minima do 
not change significantly for a given collision. 

An important point is that the interference effect de- 
scribed does not depend critically on the square barrier 
potential, and we have observed it in particular for a 
Gaussian barrier. Note that the arguments leading to 
Eq. ( |l2| ) are in fact of general validity and independent 
of the potential shape, with —d/2 and d/2 being points 
where the potential may be assumed to be essentially 
zero, and T being the corresponding transmission ampli- 
tude. The possibility to observe this effect with ultracold 
atoms rests on the ability to prepare appropriate initial 
states. Turning off the laser potential during the collision 
will leave a two peaked momentum distribution that im- 
plies at later times a visible spatial separation between 
two wave components, one faster than the other. 

We thank A. Steinberg for many useful discussions, 
and acknowledge support by Ministerio de Education y 
Ciencia (PB97-1482). 



FIG. 1. \(p\i>(t))\ 2 for differents values of t: t = (dotted- 
dashed line); t = 2.333 t u (solid line); t — 2.731t u (dashed 
line); and t = 3.233 t u (dotted line), m — 1.558023m u , Vo = 
102.5 e u , d = 2.5 l u , ~aS x = -50 l u , S x = 107.99 Z 2 , with 
an average momentum p c = 28.48 p u well above the classical 
threshold (2mV ) 1/2 = 17.87 p„. The units are scaled for 
numerical convenience in the computations as e u = 10 -13 
a.u. of energy, p u = 10~ 4 a.u. of momentum, l u = 2 x 10 6 
a.u. of lenght, m u = 10° a.u. of mass, and t u — 2 x 10 d a.u. 
of time. 

FIG. 2. Imaginary versus real parts of the incident con- 
tribution to il>iT(p,t) (empty circles), and of the transmision 
contribution (filled circles), for t = 2.731 t u and different val- 
ues of p equally spaced between p = 28 p u and p = 29 p u - 
Other parameters as in Fig. 1. 

FIG. 3. Imaginary versus real parts of u)(ui) (empty cir- 
cles) and —uj(-ut) (filled circles) for t = 2.731 t u and dif- 
ferent values of p equally spaced between the two peaks of 
Fig. 1, see the text. The prefactors corresponding to u>{ui) 
and — ut) in Vjt(P>*) are a l so shown for the same mo- 
mentum interval, solid and dashed lines respectively. Other 
parameters as in Fig. 1. 

FIG. 4. Integration contour in the complex p'-plane when 
the SDP crosses the pair of structural poles p\ and p' T . Also 
shown the structural pole p' R , the saddle point of equation 
([io|), and the incident average momentum p c . 

FIG. 5. |{p|?/>(£))| 2 as afunction ofp, for two different values 
of Vo: 102.5 e„ (solid line), and 105 e u (dashed line). The 
value of t is selected to get the maximum effect, G q ~ 0.24; 
t — 2.731 t u . Other parameters as in Fig. 1. 

FIG. 6. Imaginary versus real parts of the incident con- 
tribution to ipiriPyt) ( em Pty circles), and of the transmision 
contribution (filled circles), for Vo = 105 e„, the value of t 
for which the effect is maximum (t = 2.731 t u ), and different 
values of p. Other parameters as in Fig. 1. 
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